On the Gaussian approximation of vector-valued multiple integrals 
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Abstract: By combining the findings of two recent, seminal papers by Nualart, Peccati and 
Tudor, we get that the convergence in law of any sequence of vector-valued multiple integrals F n 
towards a centered Gaussian random vector N, with given covariance matrix C, is reduced to just 
the convergence of: (i) the fourth cumulant of each component of F n to zero; (u) the covariance 
matrix of F n to C. The aim of this paper is to understand more deeply this somewhat surprising 
phenomenom. To reach this goal, we offer two results of different nature. The first one is an 
explicit bound for d(F, N) in terms of the fourth cumulants of the components of F, when F is a 



Revalued random vector whose components are multiple integrals of possibly different orders, N 
C/3 ■ is the Gaussian counterpart of F (that is, a Gaussian centered vector sharing the same covariance 

\ with F) and d stands for the Wasserstein distance. The second one is a new expression for the 

cumulants of F as above, from which it is easy to derive yet another proof of the previously quoted 

■ result by Nualart, Peccati and Tudor. 

f-n \ Keywords: Central limit theorem; Cumulants; Malliavin calculus; Multiple integrals; Ornstein- 

■ Uhlenbeck Semigroup. 
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1 Introduction 



! Let B = (B t ) te [Q t T] be a standard Brownian motion. The following result, proved in [HE], 

yields a very surprising condition under which a sequence of vector- valued multiple integrals 
converges in law to a Gaussian random vector. (If needed, we refer the reader to section [2] 
for the exact meaning of Jj T ^ q f(t\, . . . , t q )dB tl . . . dB tq .) 

•rH , 

Theorem 1.1 (Nualart-Peccati-Tudor) Let q^, . . . , q\ ^ 1 be some fixed integers. Con- 
sider a M -valued random sequence of the form 

Fn (-^1,71) • • • j Fd t n) 

= I / fi,n(ti, ■ ■ ■ , t qi )dB tl . . . dB tqi , / fd,n(ti, ■ ■ ■ , t qd )dB tl . . . dB tqd ) , 
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where each f in G L 2 ([0,T} qi ) is supposed to be symmetric. Let N ~ ^(0,(7) be a centered 
Gaussian random vector on IR d with covariance matrix C. Assume furthermore that 

lim E[F i n Fj n ] = Cy for all i,j = l,...,d. (1.1) 

Then, as n — > oo, the following two assertions are equivalent: 

(i) F n ^ N; 

(ii) ¥/ 1 d: E[F* n ] — 3E[F? n ] 2 ->• 0. 

This theorem represents a drastic simplification with respect to the method of moments. 
The original proofs performed in [TJ [8| are both based on tools coming from Brownian 
stochastic analysis, such as the Dambis, Dub ins and Schwarz theorem. In [B], Nualart and 
Ortiz-Latorre gave an alternative proof exclusively using the basic operators 5, D and L 
of Malliavin calculus. Later on, combining Malliavin calculus with Stein's method in the 
spirit of [1], Nourdin, Peccati and Reveillac were able to associate an explicit bound to 
convergence (i) in Theorem 11.11 

Theorem 1.2 (see [4|) Consider a M. d -valued random vector of the form 
F = (F 1 ,...,F i ) 

= ( f fx(ti, ■ ■ ■ , t qi )dB tl . . . dB tqi ,...,/ f d (t u . . . , t qd )dB tl . . . dB t ) , 

\J[0,T]n J[0,T\ q d / 

where qi, . . . ,q d ^ 1 are some given integers, and each fa G L 2 ([0,T] qi ) is symmetric. Let 
C = {Cij)i^ij^a £ A^d(M) be the covariance matrix of F, i.e. = E[FiFA. Consider a 
centered Gaussian random vector N ~ ^Ki{0, C) with same covariance matrix C. Then: 

di(F,iV):= sup \E[h(F)]-E[h(N)}\ < H^ 1 || op ||C|| V 2 A C (F), (1.2) 

feGLip(l) 

with the convention ||C _1 || op = +oo whenever C is not invertible. Here: 

- Lip(l) is the set of Lipschitz functions with constant 1 (that is, the set of functions 
h : R d -»■ R so that \h(x) - h(y)\ < \\x - y\\ Rd for all x, ye R d ), 

- ||C||oj3 — su PxGK d \{o} 1 1 C* ^ 1 1 M rf / 1 1 ^ 1 1 K d denotes the operator norm on Ai d (R), 

- the quantity Ac(F) is defined as 



A C (F) :-- 



c ij ~ -(DF it DF^ m) ) 
1j / 



;i.3) 



where D indicates the Malliavin derivative operator (see section^} and (•, -)l 2 ([o,t]) is the 
usual inner product on L 2 ([0, T]). 
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When the covariance matrix C of F is not invertible (or when one is not able to check 
whether it is or not), one is forced to work with functions h that are smoother than the one 
involved in the definition (II. 2ft of di(F,N). To this end, we adopt the following simplified 
notation for functions h : E d — > E belonging to C 2 : 



\h"\ 



max sup 

i,j=l,...,d 



d 2 h 



dxidx-i 



[x 



:i-4) 



Theorem 1.3 (see [2J) Let the notation and assumptions of Theorem \1. 6 2\ prevail. Then: 

d 2 (F,N):= sup \E[h(F))-E[h(N)]\^±Ac(F), (1.5) 

||/i"|| 00 <i 1 

with Ac(F) still given by U.3\) . 

Are the upper bounds fjl.2p ~ fjl.5p in Theorems 11.21 and 11.31 relevant? Yes, very! Indeed, 
we have the following proposition. 



Proposition 1.4 (see |6|) Let the notation and assumptions of Theorem 1 1 . 1\ prevail. Re- 
call the definition lil.3\) . Then, as n — >• oo, Ac(F n ) — > if and only if E[F* n ] — 3E[F 2 n ] 2 — > 
for all i. 

In the present paper, as a first result we offer the following quantitative version of Propo- 
sition 11.41 

Theorem 1.5 Let the notation and assumptions of Theorem \1.2\ prevail, and recall the 
definition (T77J) of A C (F). Then: 

A C (F) <: j> (E[F*\ - 3E[F 2 ] 2 , E[F 2 ], . . . , E[F*] - 3E[F 2 ] 2 , E[F 2 ]) , (1.6) 

with ip : (E x E+) d E the function defined as 



d [2r\ d f 

i>(x 1 ,y 1 ,...,x d ,y d ) = ^ft-^-K 2 Y1 ( r )W V2 + ^V%N 1/4 

i,j=l \ r=l ^ ' i,j=l ^ 

+ E ^2{ qi + qj -2r)\^y^ 2 \. (1.7) 

Since, for each compact B C (0, oo) d , it is readily checked that there exists a constant 
CB, qi ,...,q d > so that 

d 

SUp ^(Si, J/i, • • -,X s ,y d ) ^ C B)?1 ,... )?d E (ki| 1/4 + W 1/2 ) ) 
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we immediately see that the upper bound (11.61) . together with Theorem II .3} now show in 
a clear manner why (ii) implies (i) in Theorem 11.11 



In a second part of this paper, we are interested in 'calculating', by means of the basic 
operators D and L of Malliavin calculus, the cumulants of any vector-valued functional 
F of the Brownian motion B. (Actually, we will even do so for functionals of any given 
isonormal Gaussian process X). In fact, this part is nothing but the multivariate extension 
of the results obtained by Nourdin and Peccati in [3J. 

Then, in the particular case where the components of F have the form of a multiple 
Wiener- ltd integral (as in Theorem II .2ft . our formula leads to a new compact representation 
for the cumulants of F (see Theorem 11.61 just below), implying in turn yet another proof 
of Theorem 11.11 (see section |4"1 



Theorem 1.6 Let m E W l \ {0} with \m\ ^ 3. Write m = l\ + . . . + h m \, where U £ 
{ei,...,e<i} for each i. (Up to possible permutations of factors, we have existence and 
uniqueness of this decomposition of m.) Consider a ~R d -valued random vector of the form 

F = (F 1 ,...,F i ) 

= ( [ fi(h,..., t qi )dB tl . . . dB tqi ,...,[ f d (t u . . . , t qd )dB tl . . . dB t ) , 

\J[0,T]n J[0,T]«d / 

where qi, . . . , ^ 1 are some given integers, and each fi £ L 2 ([0, T] Qi ) is symmetric. When 
Ik = Cj, we set Afc = j, so that F lk = F\ k for all k = 1, ... , \m\. Then: 

K m (F) = (gA w )!(H-l)! C ^( r2 ' • • ■ » r M-lK/Ai®ra/Aa ■ ■ ■ ®>Vl-i > h\ m \ ) L 2([0,T] 9A M ) ' 

where the sum ^2 runs over all collections of integers r 2 , . . . , ri m i_i such that: 
(i) 1 ^ Ti ^ q\ i for all i = 2, . . . , \m\ — 1; 

{%%) r 2 + . . . + r| m |_i — 2 7 



{%%%) r 2 < \ 2 , . . ., r 2 + . . . + r\ m \-2 < 2 ' ' . 



2 , ■ ■ ■ , ' z i • • • i ' \m\—z ^ 2 

(iv) r 3 ^ q Xl + q X2 - 2r 2 , . . ., r H „i ^ q Xl + q Xlm{ _ 2 - 2r 2 - . . . - 2r| m |_ 2 ; 
and where the combinatorial constants c q j(r 2 , . . . , r s ) are recursively defined by the relations 



and, for s ^ 3, 

c q ,i(r 2 ,...,r s ) = q Xs (r s -l)\ 



q Xl + . . . + q Xs - 2r 2 - . . . - 2r s 

r s - 1 



( ! ) 
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The rest of the paper is organized as follows. Section 2 gives (concise) background 
and notation for Malliavin calculus. The proof of Theorem 11.51 is performed in Section 
3. Finally, Section 4 is devoted to the study of cumulants, and contains in particular the 
proof of Theorem 11.61 



2 Preliminaries on Malliavin calculus 

In this section, we present the basic elements of Gaussian analysis and Malliavin calculus 
that are used throughout this paper. The reader is referred to [5] for any unexplained 
definition or result. 

Let ft be a real separable Hilbert space. For any q ^ 1, let fj® 9 be the gth tensor 
power of fj, and denote by Sj Gq the associated gth symmetric tensor power. We write 
X = {X(h), h G ij} to indicate an isonormal Gaussian process over fj (fixed once for all), 
defined on some probability space (fl, J 7 , P). This means that X is a centered Gaussian 
family, whose covariance is given by the relation E [X(h)X(g)\ = (h,g)^. We also assume 
that T = (j(X), that is, J 7 is generated by X. 

For every q ^ 1, let H, q be the qth Wiener chaos of X, defined as the closed linear 
subspace of L 2 (Q, J 7 , P) generated by the family {H q (X(h)),h G ft, ||/i|L = 1}, where H q 
is the gth Hermite polynomial given by 

We write by convention Hq = IR. For any q ^ 1, the mapping I q (h® q ) = q\H q (X(h)) can 
be extended to a linear isometry between the symmetric tensor product $) &q (equipped 
with the modified norm a/^||-||^® 9 ) and the qth Wiener chaos l-i q . For q = 0, we write 
io(c) = c, c 6 1. For q — 1, we have h(h) = X(h), h G fj. Moreover, a random variable 
of the type I q (h), h G i} 09 , has finite moments of all orders. 

In the particular case where = L 2 ([0,T]), one has that (B t ) te [ 0t j^ = (X(l[ 0it ])) te j ^ 

is a standard Brownian motion. Moreover, S) Gq = L 2 s ([0,T] q ) is the space of symmetric 
and square integrable functions on [0, T] q , and 



W)=: / f(t 1 ,...,t q )dB tl ...dB tq , fES) &q 

J[0,T]i 



coincides with the multiple Wiener-Ito integral of order q of / with respect to B, see [5] 
for further details about this point. 

It is well-known that L 2 (VL) := L 2 (n,J r ,P) can be decomposed into the infinite or- 
thogonal sum of the spaces l-i q . It follows that any square integrable random variable 
F G L 2 (Q) admits the following so-called chaotic expansion: 

oo 
<7=0 
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where f = E[F], and the f q G $) &q , q ^ 1, are uniquely determined by F. For every 
g ^ 0, we denote by J q the orthogonal projection operator on the gth Wiener chaos. In 
particular, if F G L 2 (Q) is as in (12. 8p . then J g F = I q (f q ) for every g ^ 0. 

Let {e^lfc^i be a complete orthonormal system in fj. Given / G f) 0p and g G ij 09 , for 
every r = 0, . . . ,p A g, the contraction of / and g of order r is the element of ^®(p+ < ?- 2r ) 
defined by 

oo 

f®r9= ^ • • " ^ e *r)a (8r ® (^l ® • • • ® e tr)fl® r - ( 2 - 9 ) 

U,...,«r = l 

Note that the definition of / ® r 9 does not depend on the particular choice of {ek}k^i, 
and that / £g> r g is not necessarily symmetric; we denote its symmetrization by f® r g G 
^o(p+g-2r)^ Moreover, f ®q g = f ® g equals the tensor product of / and g, whereas 
f ® q 9= (f, g)mi whenever p = q. 

It can be shown that the following product formula holds: if / G $j Qp and g G Sj Qq , 
then 



We now introduce some basic elements of the Malliavin calculus with respect to the 
isonormal Gaussian process X. Let S be the set of all cylindrical random variables of the 
form 

F = g(X(4 >1 ),...,X(4 >n )), (2.11) 

where n ^ 1, g : R n — >■ R is an infinitely differentiable function such that its partial 
derivatives have polynomial growth, and each (pi belongs to $). The Malliavin derivative of 
F with respect to X is the element of L 2 (Q,Sj) defined as 

In particular, DX(h) = h for every h G -fj. By iteration, one can define the mth derivative 
D m F, which is an element of L 2 (Q,$) @m ), for every m ^ 2. For m ^ 1 and p ^ 1, D m ' p 
denotes the closure of 5 with respect to the norm || • || mp , defined by the relation 

m 

ll^ll P m , P = SHFH + ^^OID'FIIJJ. 

i=l 

One also writes ©°° = Hm>i fl p >i ^ >m ' p - The Malliavin derivative .D obeys the following 
chain rule. If if : R n — y R is continuously differentiable with bounded partial derivatives 
and if F — (F\, . . . , F n ) is a vector of elements of D 1 ' 2 , then <p(F) G D 1,2 and 

D v{F) = Y<^ F ) DF i- ( 2 - 12 ) 



6 



The conditions imposed on (p for (I2.12p to hold (that is, the partial derivatives of ip must be 
bounded) are by no means optimal. For instance, the chain rule combined with a classical 
approximation argument leads to D(X(h) m ) = mX(h) m ~ l h for m ^ 1 and h G $). 

We denote by 5 the adjoint of the operator D, also called the divergence operator. A 
random element u G L 2 (Q,Sj) belongs to the domain of 5, noted Dom<5, if and only if it 
verifies \E(DF, u)sj\ ^ c u \\F\\l2^ for any F G D 1 ' 2 , where c u is a constant depending only 
on u. If u G Dom<5, then the random variable 8{u) is defined by the duality relationship 

E[F8{u)\ — E(DF : u)sj, (2.13) 

hl>2 



which holds for every F G 

^pc\ as J , = V 

"3 



The operator L is defined as L = YlT=o ^Q^q- The domain of L is 



DomL = {F E L 2 (Q) : ^q 2 E[(J q F) 2 } < oo} 



- jD 2 ' 2 . 

9=1 



There is an important relation between the operators D, 5 and L. A random variable F 
belongs to D 2,2 if and only if F G Dom (5D) (i.e. F G D 1,2 and DF G Dom<5) and, in this 
case, 

5DF = -LF. (2.14) 

For any F G L 2 (Q), we define L~ l F = Yl^Li ~\Jq(F)- The operator L _1 is called the 
pseudo-inverse of L. Indeed, for any F G L 2 (Q), we have that L~ l F G DomL = IB) 2 ' 2 , and 

LL- 1 F = F - E[F]. (2.15) 

We end up these preliminaries on Malliavin calculus by stating a useful lemma, that is 
going to be intensively used in the forthcoming Section HI 

Lemma 2.1 Suppose that F G D 1 ' 2 and G G L 2 (VL). Then, L^G G © 2 ' 2 and we have: 

E[FG] = E[F]E[G] + E[(DF, -DL^G)^]. (2.16) 

Proof. By fl2~Tij) and fl2T5|) . 

E[FG] - E[F]E[G] = E[F(G - E[G])] = E[F x LL~ l G] = E[F5(—DL~ l G)}, 

and the result is obtained by using the integration by parts formula f)2.13p . ■ 



3 Proof of Theorem 11.5 



The aim of this section is to prove Theorem 11.51 We restate it here for convenience, by 
reformulating it in the more general context of isonormal Gaussian process rather than 
Brownian motion. 
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Theorem 11.51 Let X = {X(h), h 6 ij} be an isonormal Gaussian process, and qd, • ■ • ,Qi ^ 

1 be some fixed integers. Consider a M. d -valued random vector of the form 

F — (Ft, . . . ,F d ) — (J ?1 (/i), . . . ,Iq d (fd)), 

where each f\ belongs to $) &qt . Let C = (CV/)i^ij<cZ -Md(R) be the covariance matrix of 
F, i.e. Cij = E[FiFj], and consider a centered Gaussian random vector N ~ Ad(0,C) 
with same covariance matrix C . Then 

A C (F) ^ V (E[Ff] " 3£[i?] a , E[F?}, . . . , E[Fj] - 3E[F 2 d } 2 , E[F 2 d }) , (3.17) 

with Aq{F) given by MJA) , and where ip : (R X — >■ R is the function given by U.7\) . 



In order to prove Theorem ll.5[ we first need to gather several results from the existing 
literature. We collect them in the following lemma. We freely use the definitions and 
notation introduced in sections [I] and [2j 

Lemma 3.1 Let F = I p (f) and G = I q (g), with f e S) Qp and g e F) Qq (p,q>l )■ 
1. If p = q, one has the estimate: 



E 



E[FG]--(DF,DG) fi 



V 



(3.18) 



r=l 



2 P" 1 



r - 1)!' 



J (2p - 2r)\(\\f ® p _ r /||J»ar + \\g ® P -r g\\l^r), 



whereas, if p < q, one has that 

2 

'Sj 

J2 P' 1 



E 



~(DF,DG) S 



^p\ 2 



q-1 
p — 1 



1)! 



2 P 



r 



{q-py-\\f\\%®v\\g®q- P g\\m 2 p 



A (p + q-2r)\(\\f® p - r f\\l®2r + 



(3.19) 



\-r S l ||j5®2r) • 



2. One has the identity: 



E[F A ] - 3E[F 



2l2 



p-1 



11/ Or /|| 

^®2p-2i- 



2p- 2r 

p — r 



||/®r/|| 

£®2p-2r 



(3.20) 



Proof. Inequalities (13 . 18[) - (13. 19[) correspond to jU Lemma 3.7] (see also [6j Proof of Lemma 
6]), whereas identity (13.201) is shown in [7J page 182]. However, for convenience of the 
reader (and also because the notation used in [7] is not exactly the same than our), we 
provide here a detailed proof of ( I3.18p . (13.1 9p and ( I3.20p . 
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1. Thanks to the multiplication formula (I2.10p . we can write 



p — 1\ f q — 1 



p/\q-± , 

pq r! ( 



r=0 
pAq 



r 



Ip+q-2-2r{f ®r+l9) 



P«E(r-l)!^_J) (J_[)w«r(/»rff)- 



It follows that 



a-^(Z)F,Z)G) fl j (3.21) 

« 2 +p 2 e: = i(^ - i)' 2 (r;)'(r!) z (p + ? - zoin/^nj.^, if P < q , 

(a - E[FG}) 2 +p 2 Zr=l(r ~ I)! 2 ( P r Z\)\ 2 P ~ 2r)!||/® r y||J aoe ^_ ap) if p = g. 
If r < p ^ g then 

ll/®r5 , |||®(p+ 9 -2r) ^ ||J Or g\\%®{v+q-2r) = (f ®p-r /, 9 ®q-r g) Sj®*r 
^ ||/®p-r/||j5»ar||^(gl 9 _ r ^ear 

^ \{\\f® P -rf\\l®2r + \\g® q - r g\\l®2r). (3.22) 

If r = p < g, then 

Wf®p9\\%^-P) < 11/ ®p »||Jb( 9 -p) < ll/IIJ®p||5'®g-pfl'llfi® 2 p- (3-23) 

By plugging these two inequalities into (I3.2ip . we deduce both (I3.18P and ( I3.19p . 

2. Without loss of generality, in the proof of ( I3.20p we can assume that fj is a L 2 -space 
of the form = L 2 (A, A, p). Let a be a permutation of {1, ... , 2p} (that is, a G © 2p ), and 
let / G $) &2p . If r G {0, . . . ,p} denotes the cardinality of {cr(l), . . . , <r(p)} D {1, . . . ,p} then 
it is readily checked that r is also the cardinality of {cx(p + 1), . . . , cr(2p)} D {p+ 1, . . . , 2p} 
and that 



/(£]., • • • , t p )f(t a (i), . . . , t(r(p))/(tp+l, • • • , t2p)f(t (T (p+l), • • • , t a (2p))dp(ti) . . . dp(t2p) 
f <S>r f(xi, • • • ,X 2p -2r) 2 dp(x 1 ) . . . dp(x 2 p-2r) = \\f ®r /|||®{2p-2r) • (3.24) 

2 



Moreover, for any fixed r G {0, . . . ,p}, there are (^) (p!) 2 permutations cr G ©2p such that 
#{cr(l), . . . , <j{p)} fl {1, . . . ,p} = r. (Indeed, such a permutation is completely determined 
by the choice of: (a) r distinct elements x±, . . . , x r of {1, . . . ,p}; (b) p — r distinct elements 
x r+ i, . . . , x p of {p + 1, . . . , 2p}; (c) a bijection between {1, . . . ,p} and {a;i, . . . , x p }; (d) a 
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bijection betwenn {p + 1, . . . , 2p} and {1, . . . , 2p} \ {x±, . . . , x p }.) Now, observe that the 
symmetrization of / ® / is given by 

f®f(ti, . . . ,t 2p ) = —— ^2 f(tcr{i),---,ta{p))f(ta( P +t),---,t cr (2p))- 

Therefore, 

||/®/||$®2j> = / x, 2 ^ / f(t er {l),'",ta(p))f(U{p+l),---,t<T(2p)) 

{2p)l *yee 2p jA2p 

X /(V(1), • • • ,ta'{p))f(ta'{p+l), ■ ■ ■ ,t a >(2p))dn(tl) . . . dfj,(t 2p ) 

= 7o~Ti 5^ / f(tl,---,tp)f(t p+ i,...,t 2 p) 

x f(t<j(i), ■ ■ ■ )*G-(p))/(*o-(p+i), • • • ,t a (2p))d[i{ti) . . . dfi(t 2p ) 
1 P f 

/ f(ti,---,tp)f(t p+1 ,...,t 2 p) 

»}n{l,...,p}=r 

X f(ta(l)i ■ ■ ■ ,to{p))f(ta{p+l), ■ ■ ■ ,t<r{2p))d^{t\) ■ ■ ■ dfl(t 2p ). 



(2p)\ 

{<7(l),...,<T(p)}n{l,...,p}=r 



P-1 / \ 2 



Using ( I3.24p . we deduce that 

(2p)!||/®/|li«- = 2(p!) 2 ||/||^ + (p!) 2 £ (") 11/ ®r /!£.(*-*>■ (3.25) 
On the other hand, we infer from the product formula f)2.10p that 

^ = W) 2 = X>! (f) hp-2r{f®rf). 
r=0 ^ ' 

Using the orthogonality and isometry properties of the integrals I p , this yields 

P / \ 4 



E[F 4 ] = J>!) 2 (^) (2p-2r)!||/® r /||^ (2p _ 2r) 

r=0 ^ ' 



(2p)!||/®/||J« (a p> + (p!) 2 ||/||^ + £(r!) a ( P Y(2p - 2r)!||/® r /||J 9 ^_ ar) . 



r=l 



By inserting (I3.25P in the previous identity (and because (p!) 2 11/11^®? = E[F 2 ] 2 ), we get 

dSSOD- 

□ 
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We are now ready to prove Theorem 11.51 If Z G L 4 (f2), as usual we write Xi{Z) = 
E[Z 4 } - 3E[Z 2 } 2 for the fourth cumulant of Z. We deduce from (l372Tjj) that, for all p ^ 1, 
/ G i} 0p and r G {1, . . . ,p — 1}, one has X4,(Ip(f)) ^ and 



11/ ®r /Ilfl®2p-2r ^ 



r! 2 (p — r)! 2 



P 



|4 



X4(/ P (/))- 



Therefore, if /, g G .f) 0p , inequality (I3.18P leads to 

2" 



J5[/p(/)W]-^<^W).^)>fl 



[X4(w))+x 4 (/^))]x: ;^_ 2 r ^ 

^[X4(W))+X4(W)]E 



r=l 



(3.26) 



On the other hand, if p < q, f G i} 0p and g G inequality (I3.19P leads to 

2 



E 



V 



(DI p (f),DI q (g)) 



to 










( 









.q 



(DI p {f),DI q (g)) s 



< £[W) 2 ] v^Q?)) + i E r2 o° + 9 - 2r ) ! 

r=l 



,! 2 



X4(/p(/)) 



P 



1 2 



X4(/«j(c/)) 



(g — r)! 2 p! 2 /X ' 1V pw yy ' (p — r)! 2 g! 2 

< E[i P (J)*\y/x*{i q {g)) + \ X> + ? - 2r ) ! 



X 



X4(/p(/)) 



X4(i g (»)) 



so that, if p 7^ g, / G ij 0p and g G -fj 09 , one has that both E 



l(DI p (f),DI q (g)) s 



and _E 



UDI p (f),DI q (g)) s 



are less or equal than 



p/\q—l 



(3.27) 



.-£(p + g-2r)! * 

r=l ^ ' 



X4(/p(/)) 



Since two multiple integrals of different orders are orthogonal, on has that 
Qj = E[EiEj] = E[I qi {fi)I qj {fj)} = whenever Qi ^ q r 

Thus, by using (I3.26p -( 13.27p together with sjx\ + . . . + x n Si v /x7+. . .+y/x^, we eventually 
get the desired conclusion f)3.17p . □ 
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4 Cumulants for random vectors on the Wiener space 



In all this part of the paper, we let the notation of section |2] prevail. In particular, X = 
{X(h), h G $)} denotes a given isonormal Gaussian process. 

4.1 Abstract statement 

In this section, by means of the basic operators D and L, we calculate the cumulants of 
any vector-valued functional F of a given isonormal Gaussian process X. 

First, let us recall the standard multi- index notation. A multi-index is a vector m = 
(mi, . . . , rrid) of N d . We write 

d r\ d 

i=l 1 i=l 

By convention, we have 0° = 1. Also, note that \x m \ = y m , where yi = \xi\ for all i. If 
s G N d , we say that s ^ m if and only if s, ^ rrii for all i. For any i = 1, . . . , d, we let 
Ci G M d be the multi-index defined by (ej)j = 5ij, with 5{j the Kronecker symbol. 

Definition 4.1 Let F = (Fi, . . . , Fd) be a R. d -valued random vector such that E\F\ m < oo 
for some m £ W l \ {0}, and let 4>F{t) = E[e l ^' F '» d ] , t G W 1 , stand for the characteristic 
function of F. The cumulant of order m of F is (well) defined by 

K m {F) = {-i)\ m \d m \o g Mt)\t=o- 

For instance, if F h Fj G L 2 (fi), then k £i (F) = E[Fj\ and K ei+e .{F) = Cov[Fj, Fj\. 
Now, we need to (recursively) introduce some further notation: 

Definition 4.2 Let F = (Fx, . . . , Fd) be a R. d -valued random vector with Fi G D 1 ' 2 for each 
i. Let li, / 2 , • • • be a sequence taking values in {ei, . . . , e^}. We set r^(F) = F l1 . If the 
random variable r^ lj i k (F) is a well-defined element of L 2 (Q) for some k ^ 1, we set 

T h! ..., lk+1 (F) = (DF^^-DL-'T^F))^. 

Since the square-integrability of F; lv .. i / fc (F) implies that L~ 1 r/ lv .. i / ( ,(F) G DomL C D 1 ' 2 , 
the definition of r^...^ +1 (F) makes sense. 

The next lemma, whose proof is left to the reader because it is an immediate extension 
of Lemma 4.2 in [3] to the multivariate case, gives sufficient conditions on F ensuring that 
the random variable i k {F) is a well-defined element of L 2 (Q). 

Lemma 4.3 1. Fix an integer j ^ 1, and assume that F = (Fi,...,Fd) is such that 
Fi G W' 23 for all i. Let li, Z 2 , • • • , lj be a sequence taking values in {e\, . . . , e^}. Then, 
for all k = 1, . . . , j, we have that r^ lv . i k (F) is a well-defined element of W~ k+1,23 +1 ; in 
particular, one has that Y^ j.^F) G D 1,2 C L 2 (Q) and that the quantity E^r^^^F)] is 
well-defined and finite. 

2. Assume that F = (Fi, . . . , Fd) is such that Fi G D°° for all i. Let li, I2, ■ ■ ■ be a 
sequence taking values in {ex, . . . , Cd}- Then, for all k ^ 1, the random variable Ti lr ^j k (F) 
is a well-defined element 0/D 00 . 
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We are now ready to state and prove the main result of this section, which is nothing 
but the multivariate extension of Theorem 4.3 in [3]. 

Theorem 4.4 Let m G N d \ {0}. Write m = l\ + . . . + l\ m \ where U G {ei,...,erf} 
for each i. (Up to possible permutations of factors, we have existence and uniqueness of 
this decomposition of m.) Suppose that the random vector F = (F\, . . . ,Fd) is such that 
Fi G ©l m l' 2 ' m1 for all i. Then, we have 

K m (F) = (\m\-l)\E[T h _ ll jF)]. (4.28) 

Remark 4.5 A careful inspection of the forthcoming proof of Theorem 14.41 shows that the 
quantity E^Ti limm h m (i* 1 )] in (I4.28P is actually symmetric with respect to l\, . . . , l\ m \, that is, 

Va G 6 W , E[T h _ llm] (F)] = E[r la(1) ,... M]m]) (F)]. 

Proof of Theorem \4-4\ The proof is by induction on \m\. The case \m\ = 1 is clear because 
K ej {F) = E[Fj] = E[r ej (F)] for all j. Now, assume that (1428]) holds for all multi-indices 
m G N d such that \m\ ^ N, for some N ^ 1 fixed, and let us prove that it continues to 
hold for all the multi-indices m verifying \m\ = N + 1. Let m G N d be such that \m\ ^ N, 
and fix j = 1, . . . , d. By applying repeatidely (I2.16P and then the chain rule (I2.12p . we can 
write 

E[F m+e i] = E[F m x F e .(F)] 

= E[F m ]E[Te.(F)] + E[(DF m , —DL~ l T e .(F))^] 

= E[F m ]E[T ej (F)]+ E[F m - l ^{DF l ^,-DL- 1 r e .(F)) Si } 
= E[F m }E[F e .(F)} + E[F™- l nr ej , ln (F)} 

= E[F m ]E[T ej (F)]+ J2 E[F m - l n]E[T ej!hi (F)} + E E[F m - l n- l ^T ej:l 

l^ll^|m| l^i 1 ,i 2 ^|m| 

il ,12 different 

= E[F m ]E[T ej (F)]+ E[F m - l n]E[T e]Mi (F)] 
+ E[F m - l n-^]E[r ej , lii!li2 (F)} 

,%2 tfi\ 
i\,%2 different 

l<il,...,i| m |_ 1 <|m| 
2l,...,2| m |_l pairwise different 

+|m|!^[r ej ,, ij ...,. m| (F)] 
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so that, using the induction property, 

E [F m ^\ = E[F m dn ej (F)+ £ E{F m ~ l ^K ej+hi (F) 



0! 

E ^ 



^2 



2j ~t~^»2 ^^"^ 



il different 



+ ... 



E ^ 



m l il ... Jj| m ,_ 1 



1 



l<il,...,i| m |_ 1 <|m| 
n,.--^|m|-l pairwisc different 



] (m-l)! Ke ^ + - +Z Vl-i (F) 



+|m|!E[r ei ,, i ,..., V| (F)] 
E ^[^ m "1^T«e, +s (F) #5 S + Iml^^^,...,^ (F)]. 



|s|^m— 1 



Here, _B S stands for the set of pairwise different indices ii, . . . , iy G {1, . . . , |m|} such that 
li x + . . . + /», , = s, whereas denotes the cardinality of B s . Also, let Dj = {i = 

l,...,\m\ : /j = ej} and observe that m = (mi , . . . , 777^) with = #Dj. For any 
s «C m, it is readily checked that #B S = (™ : ) . . . (Indeed, to build a multi- index 

s = (s\, . . . , Sd) so that s ^ m, one must choose s\ indices among the mi indices of D\ 
up to Sd indices among the rrid indices of Dd, and then the order of the factors in the sum 
l ix + . . . + L, .) Therefore, 



E[F 



E 



|s|^m— 1 



17li 
Si 



E[F m ~ s ] K ej+S (F) + \m\\E[T ej:lii ,..., liM (F)} 



E 



mi 

Si 



m d 

Sd 



(-z) |mh|s| 9 m " s F (O) x (-z) |s|+1 ^ +s log0 F (O) 



+|m|!E[r e . Ai ,... V| (F)]-K e . +m (F) 



i)M+ifln»(^ io g ^)( ) + ImllE^^, , (F)\ - K e . +m (F) 



dt j 

= H) |m|+1 d m +^<M0) + \m\\E[T ehki _ k ^ (F)] - K ej+m {F) 
= E[F m+ej ] + \m\\E[T ejMi _ him[ (F)) - ^ +m (F), 

leading to 

|m|!S[r e . Aii ... A|m| (F)] = K ej+m , 

implying in turn that (14.281) holds with m replaced by m + ej. The proof by induction is 
concluded. □ 
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4.2 The case of vector-valued multiple integrals 

We now focus on the calculation of cumulants associated to random vectors whose com- 
ponent are in a given chaos. In (I4.29|) (and in its proof as well), we use the following 
convention. For simplicity, we drop the brackets in the writing of f\ l ® r2 ■ ■ ■ ®r {ml _ 1 f\ lml _ 1 , 
by implicitely assuming that this quantity is defined iteratively from the left to the right. 
For instance, f® a g®ph(g)~ / k actually means ((f<^ a g)<S>/3h)<^jk. 

For convenience, we restate Theorem 11.61 (in the more general context of isonormal 
Gaussian process). 

Theorem 4.6 Let m E N d \ {0} with \m\ ^ 3. Write m = h + . . . + l\ m \, where U G 
{ei, . . . ,ed} for each i. (Up to possible permutations of factors, we have existence and 
uniqueness of this decomposition of m.) Consider a M. d -valued random vector of the form 

F = (Fi, ...,F d )= (l qi (fl), Iq d {fdj), 

where each /, belongs to $) &qi . When Ik = ej, we set Afc = j, so that F lk = F Xk for all 
k = 1, . . . , \m\. Then: 

K m {F) = (?a m )!(H - 1)! J^c^(r 2 , . . . ,r\ mhl )(f Xl ® r2 f\ 2 • • • ^r^.Jx^, fx lml ) ^x ]m] , 

(4.29) 

where the sum ^2 runs over all collections of integers r 2 , . . . , J"| m |_i such that: 
(i) 1 ^ Ti ^ q Xi for all i = 2, . . . , \m\ — 1; 
(n) r 2 + ...+ r| m |_i = lj l_l ; 

(m) r 2 < \ 2 , . . ., r 2 + . . . + r\ m \-2 < 2 ' ' ; 

(iv) r 3 ^ q Xl + q X2 - 2r 2 , . . ., r\ m \^ < q Xl + <?a m _ 2 - 2r 2 - . . . - 2r H _ 2 / 
and where the combinatorial constants c 9j ;(r 2 , . . . , r s ) are recursively defined by the relations 

c,M) = ^(r 2 -l)'-(^:l)(\ 
and, for s ^ 3, 

c>,.,r s ) = ^(r.-l)l(^ + - + ^- r ^-"- 2r - 1 - 1 



r 2 - 1 
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Proof. If / G Sj Qp and g G i} 09 (p, q ^ 1), the multiplication formula yields 
(Dlp(f), —DL^ 1 I q (g)) S j = p(I p „ x {f),I q ^(g)). 



r=0 

pAg 

-2r(/® r y). (4.30) 



Thanks to (|4.30j) . it is straightforward to prove by induction on \m\ that 

r h ,...,i lmi (F) (4.31) 

3A 1 A?A 2 [9A 1 +-+3A |m| _ 1 -2r 2 -...-2r-| m |_ 1 ]A 9A|m| 



■■■ Y c 3i i(r 2 ,...,r H )l 



r 2 =l r M =1 
{r 2 +...+r| m |_ 1 < ! } 



{^2< \ 2 } 



Xl r .. , , ,^i+- +g A| TO |_ 1 l -^A 1 +-+gA |m -2r 2 -...-2r| m | (f\i®r 2 f\2 ■ ■ ■ ®r\ m \f\\ m[ ) ■ 



(4.32) 

Now, let us take the expectation on both sides of (I4.32p . We get 
= (|m|-l)!J5[T Ili ... iI|m| (F)] 

qx,Aqx 2 [9A 1 +...+9A| m |_ 1 -2r 2 -...-2r| m |_ 1 ]A 9A|m| 

= (H-l)! ■■■ c 9 , z (^,...,r H )l {r2<!Ai +^ } ... 

r 2 =l r \m\ =1 

x ^ r2+ ... + ,„ ll - 1 < " l+ ■■r | '" | - l } 1 ^ 2+ ■■■ + ^r 9Al+ ^^r v| } x 

Observe that, if 2r 2 + . . . + 2r| m | = q Xl + • • • + <7a w and r| m | < q Xl + . . . + <2A| mM - 2r 2 - 
. . . - 2r| m |_i, then 

2r H = 9A W + (?Ax + • • • + gA| m ,_! - 2r 2 - . . . - 2r M _!) ^ g A|m| + r H , 
that is, r| m | ^ g A|m| , so that r\ m \ = <?a m • Therefore, 

qx ± Mx 2 [9A 1 +-+9A| m |_ 1 -2r 2 -...-2r-| m |_ 1 ]A 9 A |m| 

= (H- 1 ) 1 ••• Yl Cg,»(^,---,nm|)l [ri ^A 1 +«A 3 ] ••■ 



''2=1 r M =1 
X l 9A l + "' + 9A |m|-l ^ 9A l + "' + ,A lm 

{r 2 +...+r ]m |_ 1 < J - } { r 2+...+r-| m | = 5 



X (f\i®r 2 f\2 ■ ■ ■ ®r| m |_i/\| m |_u /A| m |)^®9A| m| , 

which is the announced result, since c gjZ (r 2 , . . . , r| m |_i, q\^) = (q\ M ) ] -c q ,i(r2, ■ ■ ■ ,r\ m \-i). n 
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4.3 Yet another proof of Theorem 11.11 

As a corollary of Theorem 14. 6[ we can now perform yet another proof of the implication 
(ii) —> (i) (the only one which is difficult) in Theorem 11.11 So, let the notation and 
assumptions of this theorem prevail, suppose that (ii) is in order, and let us prove that 
(i) holds. Applying the method of moments/cumulants, we are left to prove that the 
cumulants of F n verify, for all m G N d , 

, r > f if \m\ 7^ 2 

K m (F n ) -> K m (N) = | ^ if m = ei + e - aS U ~* 00 ' 

Let m G N d \ {0}. If m = for some j (that is, if and only if \m\ = 1), we have 
Km(F n ) = E[Fj )Tl ] = 0. If m = Ci + for some i,j (that is, if and only if \m\ = 2), we have 
K m (F n ) = E[Fi tTl Fj jn ] — >■ Cjj by assumption (II. ip . If \m\ ^ 3, we consider the expression 
(I4.29p . Thanks to ( I3.20p . from (ii) we deduce that ||/ i)Tl ® r /i, n ||i 2 ([o,T]'i) — ^ as n — > oo 
for all i, whereas, thanks to (jl.l[), we deduce that ||/i,Ti|||,aQo,Ti*») 

= E[F? n ] for all 

2, so that sup,^ ||/j,n||L 2 ([o.r]9i) < oo for all z. Let r% . . . ,ru_i be some integers such that 
(i)-(iv) in Theorem 14.61 are satisfied. In particular, r 2 < 9Xl „ 9X2 . From fl3.22p - fl3.23p . it 
comes that \\f\ 1 ,n®r 2 f\2,n\\ L 2^ T iix 1 +"x 2 ~ 2r2 ) — ► as n — > oo. Hence, using Cauchy-Schwarz 
inequality successively through 

\\g® r h\\ L 2([Q T ]p+ q -2r} ^ \\g ® r /l|| L 2([ 0]T ]p+ 9 -2r) ^ 1 1 # 1 1 L 2 ( [0,T]P ) | N I L 2 ( [0,T] <2 ) 

whenever g G L 2 s ([0,T] p ), h G L^O.T]") and r = 1, . . . ,p A g, we get that 

(/\ ll n®ra/A 3 ,n ■ ■ ■ ®r w .,/vi- 1 ,n> A|m|^) La( [ ,r| B Vl) ^ as ™ ^ °°- 

Therefore, K m (F n ) — >• as n — >■ oo by fl4.29p . j-j 
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